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1. Introduction 


There are many design requirements associated with tilt-rotor performance including a low disc 
loading in the hover configuration and the ability to rotate the rotors forward to achieve cruise 
speeds up to 450 knots [1]. The problem becomes more complex since in vertical flight and in 
hover, a large portion of the rotor is directly over the wing which produces a large downwash 
effect upon the wing. The downwash effect increases thrust requirement of the aircraft by 
approximately 10-12 percent [2]. Other problems associated with this configuration are related to 
high helical tip Mach numbers (M t i p ) which represent a critical performance issue in high speed 
cruise (350 - 450 knots). Aeroelastic stability is another important consideration in the design of 
tilt-rotors. Due to the large thrust requirement in hover, the prop-rotors have a much greater radius 
than standard propellers. This increases the tip speed which in cruise may cause individual blade 
flutter or a coupled flexible motion between the rotor, wing and pylon known as whirl flutter. 
Also, blade sweep which is introduced for improvements in aerodynamic performance can be 
detrimental for aeroelastic stability. Since civil tiltrotors are required to be stable at a 20 percent 
margin above their dive speed, this means that the flutter speed must be in the order of 620 knots 
for a tiltrotor with a cruise speed of 450 knots. 

There are several different techniques which can be used to address these issues. For example, 
the tip Mach number can be reduced through rotor tip speed reduction or through the use of blade 
sweep which reduces the effective Mach number. Another alternative is to increase the drag 
divergence Mach number (M^d) at the tip to value above M t j p . This can be accomplished through 
reductions in the blade thickness. However, each of these options will adversely affect the hover 
performance, drive system weight or aeroelastic stability of the rotor blade. In the helicopter 
mode, to maintain a high figure of merit in hover, the solidity of the blade must be increased since 
thinner airfoils are used for maintaining efficiency in cruise. Due to such conflicting requirements 
between hover, conversion and airplane mode prop-rotor performance, the use of formal numerical 
optimization techniques is appropriate for studying the design trade-offs associated with the 
development of high speed tilt-rotors. 



Figure 3 XV- 15 in cruise mode 
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2. Objectives 


The objective of this study is to develop an optimization procedure for high-speed and civil tilt- 
rotors by coupling all of the necessary disciplines within a closed-loop optimization procedure. 
Both simplified and comprehensive analysis codes are used for the aerodynamic analyses. The 
structural properties are calculated using in-house developed algorithms for both isotropic and 
composite box beam sections. The major objectives of this study are as follows: 

1 . Aerodynamic Optimization : Investigate the effects of blade aerodynamic characteristics on 
cruise and hover performance of prop-rotor aircraft using the classical blade element 
momentum approach with corrections for the high lift capability of rotors/propellers. 

2. Coupled Aerodvnamic/Structures Optimization : Develop a multilevel hybrid optimization 
technique for the design of prop-rotor aircraft. The design problem is decomposed into a 
level for improved aerodynamics with continuous design variables and a level with discrete 
variables to investigate composite tailoring. The aerodynamic analysis is based on that 
developed in the first objective and the structural analysis is performed using an in-house 
code which models a composite box beam. The results are compared to both a reference 
rotor and the optimum rotor found in the purely aerodynamic formulation (objective 1). 

3 . Multipoint Optimization : Extend the multilevel optimization procedure of objective 2 to a 
multipoint design problem. Include hover, cruise and take-off as the three flight conditions 
to be simultaneously maximized. 

4. Coupled Rotor/Wing Optimization : Using the comprehensive rotary wing code, 
CAMRAD, develop an optimization procedure for the coupled rotor/wing performance in 
high speed tilt-rotor aircraft. The developed procedure contains design variables which 
define the rotor and wing planforms. 

These objectives have all been achieved and the details of these investigations are presented in 
the next four sections. 
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3. Aerodynamic Optimization 


The objective of this research effort is to explore the use of a simpler computational tool limited 
to aerodynamic performance only as one level of a multilevel based optimum design problem. 
Ultimately the results of this analysis, in a complete formulation, would represent an input to a 
more comprehensive structural and aeroelastic representation of the rotor characteristics. The 
results of this research effort represent evaluation of the suitability of optimization procedures to 
improve the aerodynamic performance of high speed prop-rotors by including all aerodynamic 
parameters as design variables. For example, blade chord, twist, thickness and camber 
distributions are allowed to vary during the design optimization procedure. Since the design 
involves multiple objective functions, a multiobjective function formulation technique is used. 
Necessary side constraints are imposed on the design variables, during optimization, to prevent 
unrealistic values. However, in a multidisciplinary coupled problem, based on multilevel 
decomposition formulation, more realistic constraints on chord, twist, thickness and camber would 
be obtained from the optimal sensitivity parameters or the coordination constraint equations linking 
this formulation to the structural and aeroelastic representations in the other optimization levels. 

3.1 Optimization Problem 

The objectives of this study are to maximize the propulsive efficiency in high speed cruise (T| c ) 
and the figure of merit in hover (FM). Two different forward speeds are used to represent the high 
speed cruise flight condition corresponding to 300 knots and 400 knots. In case of the 300 knots 
forward speed condition, the optimization is first performed individually for both hover and cruise. 
Then a formal multiobjective function formulation technique is used to simultaneously maximize 
the hover figure of merit and cruise propulsive efficiency. The multiobjective formulation 
technique is also used in case of the 400 knots forward speed. Constraints are imposed on the 
rotor thrust in both hover and in high speed cruise. The blade is discretized and the values of 
chord (c), twist (9), thickness to chord ratio (t/c) and zero lift angle of attack (otzl) at each segment 
are used as design variables. The optimization problem is summarized as follows. 
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Maximize 


FM(O), Tic(O) 


subject to Th = (Th)ref 

T c = (T c ) re f 

where O = [c(y), 6(y), a z i(y), t/c(y)] is the design variable vector and the subscript ‘ref is used to 
indicate reference or baseline rotor values. 

3.2 Multiobjective Formulation 

Since the optimization problem involves more than one design objective, the objective function 
formulation is more complicated. In most of the existing work, the individual objective functions 
are combined using weight factors in a linear fashion. Such methods are judgmental as the answer 
depends upon the weight factors which are often hard to justify. Therefore, the problem is 
formulated using the Kreisselmeier-Steinhauser (K-S) function approach [3]. Using this function 
the multiple objective functions and constraints are transformed into a single envelope function 
which is then minimized using unconstrained optimization techniques. The K-S function has been 
found to perform extremely well by the authors for multiobjective rotary wing optimization 
problems [4,5]. 

The first step in formulating the objective function in the Kreisselmeier-Steinhauser (K-S) 
function approach involves transformation of the original objective functions into reduced objective 
functions. When the individual objective functions are to be minimized, these reduced objective 
functions assume the following form. 

Flcs) = - l.o - g m » £ o k = 1 NF min (la) 

When the individual objective functions are to be maximized, the reduced objective functions are as 
follows. 

FJJ(O) = 1.0 - - gmax ^0 k = 1, ... , NF max ( lb ) 
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where Fk 0 represents the value of the original objective function Fk calculated at the beginning of 
each iteration and O is the design variable vector. The quantity g max is the value of the largest 
constraint corresponding to the original constraint vector, gj(<l>) (j = 1, 2, ••• , NC), and is held 
constant during each iteration. These reduced objective functions are analogous to constraints, 
therefore a new constraint vector f m (0) (m = 1, 2, — , M where M = NC + NF) is introduced 
which includes the original constraints and the constraints introduced by the reduced objective 
functions (Eqns. 1). The design variable vector remains unchanged. The new objective function 
to be minimized is defined using the K-S function as follows. 

1 M 

F KS(<J>) = fmax + g In ^exp[5(f m - f max )] (2) 

m=l 

where f max is the largest constraint corresponding to the new constraint vector, f m (0), and in 
general is not equal to g max . The optimization procedure is as follows. Initially in an infeasible 
design space, where the original constraints are violated, the constraints due to the reduced 
objective functions (Eqns. 1) are satisfied, i.e. g max is negative. Once the original constraints are 
satisfied, the constraints due to the reduced objective functions become violated (positive). When 
this happens, the optimizer attempts to satisfy these constraints introduced by the reduced objective 
functions and in doing so, minimizes the original objective functions (F^). The multiplier 8 is 
analogous to a draw-down factor where 5 controls the distance from the surface of the K-S 
objective function to the surface of the maximum constraint function. When 8 is large, the K-S 
function closely follows the surface of the largest constraint function. When 8 is small, the K-S 
function includes contributions from all violated constraints. Although in this traditional 
formulation, the objective functions and constraints are coupled in a less judgmental manner, it 
must be noted that weight factors can be introduced for the objective functions and/or constraints 
through Eqns. 1. For a given set of weighting factors, the K-S function only selects one point on 
a locus of optimum points on a plot of hover figure of merit versus propulsive efficiency. 
However, as in linear combinations of objective functions, weight factors can be used in the K-S 
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formulation to emphasize the importance of one or more objective functions over the remaining 
objective functions. In the present formulation, no such weight factors are used. 

3.4 Analysis 

The aerodynamic formulation is based on the two dimensional compressible aerodynamic 
representation developed by Smith [6] and later used by Talbot [7] in the formulation of an axial 
flow performance analysis. The modifications by Talbot include an empirical correction to two 
dimensional stall behavior in order to represent the high lift capability demonstrated by rotors and 
propellers. Compressibility and Reynold’s number effects are included in this formulation. The 
blade element theory used in the algorithm is due to Glauert [8]. In this formulation, the sectional 
lift and drag are resolved into elemental thrust and torque for each section of the blade. The force 
and momentum equations for thrust and torque assume the following form. 

dTi = 47trpoo(Voo + Di)\)idr 

dT? = -^pooW 2 c[cicos(ai) - Cdsin(ai)]dr 
2 

dQi = 47tr 2 pc»(Voo + vOujdr 

dQ 2 = -pooW2c[cisin(ai) + cdcos(oci)]rdr 

2 

where dT, dQ and dr represent the section thrust, torque and element length, respectively, W is the 
resultant velocity, V<» is the forward velocity and 'Oi and uj represent the inflow and swirl 
velocities, respectively. The chord length and radial locations are denoted c and r, respectively and 
Poo is the density of air. The sectional coefficients of lift and drag are C| and Cd, respectively and oq 
is the local inflow angle of the blade section. The subscripts (1) and (2) correspond to the 
momentum and force equations, respectively, for thrust and torque. Detailed expressions for the 
coefficients of lift and drag (cj and Cd), which represent the high lift capability of rotary wings in 
post stall angle of attack region, are found in Ref. 6. This system of equations is then used to 
solve for the inflow and swirl velocities by equating the thrust and torque as follows: 


(3) 

(4) 

(5) 

( 6 ) 



In the initial study performed by Smith [6] an empirical fit was performed on NACA 63 and 64 
series airfoil families in order to supply a functional relationship between maximum lift coefficient 
and the sectional thickness and camber for incompressible flow. These functional relationships 
were later modified by Talbot [7] to model the Advanced Tiltrotor Blades (ATB) [9,10]. This 
representation is suitable only for conceptual design and it is recognized by the authors that a more 
sophisticated airfoil property representation is desirable during detailed design. However, as the 
baseline XV- 15 rotor system utilizes airfoils of this general capability, it is felt that this 
representation is adequate to demonstrate the capability of an optimization procedure to generate a 
feasible design. The analysis is easily implemented within an optimization procedure and offers 
significant computational advantages compared to more comprehensive analysis procedures. The 
adequacy of this representation is demonstrated in Figs. 4 and 5, where this relatively simplistic 
formulation is shown to correlate very well with measured axial flow performance of the XV- 15 
rotor system in both hover and in airplane mode from tests conducted at the Outdoor Aerodynamic 
Research Facility (O.A.R.F.) and from wind tunnel test data obtained at NASA Ames Research 
Center [9,10]. This representation of the rotor, which is representative of the original design point 
of the XV- 15 tiltrotor, is used as the baseline, or reference, rotor for both the 300 and 400 knots 
cases. Further, the results obtained using this approach are comparable with results obtained by 
Dadone et al. using comprehensive rotorcraft analyses [1]. Other elementary algorithms have been 
proposed to model the post stall delay due to rotation [11]. 

3.5 Optimization Implementation 

A nonlinear programming procedure, based on the method of feasible directions, is used for 
the optimization [12]. During optimization, several evaluations of the objective function and 
constraints are necessary to obtain a feasible design. Therefore, the process can become 
computationally expensive if exact analyses are performed for every function evaluation. An 
approximate analysis technique is therefore implemented in the calculations of the objective 
functions and the constraints. The approximate analysis used in this study is based on a two point 
exponential procedure [13]. 
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3.6 Results and Discussions 


The reference rotor used is the original XV- 15 prop-rotor [9,10]. In both airplane mode 
forward flight conditions the rotational speed is 421 RPM and the aircraft altitude is 25,000 feet. A 
vehicle weight of 13,000 lb. and aircraft lift to drag ratio (L/D) of 8.4 is assumed for both cruise 
flight conditions. Therefore, the thrust in cruise is constrained to be 774 lb. for the two engine 
aircraft. In hover, the aircraft is assumed to be operating at sea level conditions with a rotational 
speed of 570 RPM and a 12 percent down load effect from the rotor/wing interaction. The thrust 
in hover for the two rotor aircraft is therefore constrained to be 7280 lb. To ascertain the credibility 
of the optimization procedure as a design tool, the first problems addressed are the optimization of 
the XV- 15 prop-rotor, initially for hover only and for 300 knots cruise operation, next. The intent 
of these single objective function optimization procedures is to demonstrate the capability of the 
optimization procedure in identifying the planform characteristics known to be optimum for both a 
helicopter in hover and a propeller. The multiobjective function formulation optimization 
procedure is then performed first for the 300 knots cruise condition and then for the 400 knots 
forward flight condition. These flight conditions are summarized in Table 1. The blade is 
discretized into 10 segments. The discrete nodal values of chord (c), twist, relative to the twist at 
the y = 0.75 radial location, (0), thickness to chord ratio (t/c) and zero lift angle of attack (a z i) are 
used as design variables yielding a total of 44 design variables. The following side constraints are 
imposed on the design variables throughout the blade span. 


0.001 

< 

c/R 

< 

0.50 

(9a) 

-15° 

< 

0 

< 

45.0° 

(9b) 

0.05 

< 

t/c 

< 

0.40 

(9c) 

-5.0° 

< 

Otzl 

< 

3.0° 

(9d) 


It must be noted that the lower bound of c/R = 0.001 listed above (9a) is included to consider the 
aerodynamic only optimization case. Although the same values are used as side constraints for the 
chord in the other optimization cases as well so that meaningful comparisons of the chord 
distributions can be made, the actual chord values in these cases never approach this bound. 
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Table 1 Summary of Flight Conditions (Objective 1) 


Vehicle Weight 

13,000 lb 

Blade Radius 

12.5 ft 

Hover Altitude 

Sea level 

Hover Thrust, Th 

7280 lb 

Hover Rotational Speed 

570 RPM 

Cruise Altitude 

25,000 ft 

Cruise Thrust, T c 

7741b 

Cruise Rotational Speed 

421 RPM 


The optimization results for the following cases: hover, airplane mode with forward speed of 
300 knots and multiobjective formulation (hover and cruise), also at 300 knots, are presented in 
Figs. 6-10. Figure 6 shows that the single objective function formulation produces designs with 
significant performance improvements. In case of hover, the figure of merit is increased by 9.2 
percent and in case of airplane mode the propulsive efficiency is increased by 8.7 percent. In the 
multiobjective function case, as expected, the improvements are significant, but not as dramatic as 
in the single objective function cases. In this case, the Figure of merit is increased by 7.7 percent 
and the propulsive efficiency is increased by 3.7 percent. It is important to note that the 
multiobjective formulation is capable of retaining the bulk of the hover performance improvement, 
but the propulsive efficiency is apparently limited by the solidity requirements for hover. 

Examinations of the basic planform generated by this procedure in Fig. 7 shows that the 
optimization procedure is capable of identifying the classic optimal planform for both hover and 
cruise operations. In the case of hover, the planform varies as the inverse of the radial location and 
in the cruise case (propeller mode), the root chord is driven almost to zero. As shown in the 
figure, the multiobjective function planform is more closely related to the hover planform although 
the airplane mode efficiency objective function does reduce the chord values from root to midspan. 
The reason the multiobjective planform more closely models the hover only solution is due to the 
fact that a much larger solidity is necessary to achieve the required thrust in hover. The outboard 
chord reduction, relative to the reference, in the multiobjective design is a reflection of the lack of a 
maneuver margin requirement in the optimization problem formulation. Figure 7 not only shows 
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the optimizer’s ability to recognize the classical solution, but also the root chord values in the case 
of the cruise only optimization demonstrate the necessity of imposing realistic constraints on the 
design variables during the optimization process. It is obvious more practical side constraints are 
necessary for a realistic design. 

□ Reference 

Single Objective Function 
9 Multiobjective Function 


8 . 7 % 



Figure 6 Summary of optimum results 


The twist (0) distributions for the reference and the three optimum rotors are presented in 
Fig. 8. In case of the hover only optimization, the classical solution, which varies as the inverse 
of the radial location, is once again achieved. In the airplane mode, the twist distribution is much 
more linear which is expected due to the high inflow velocity at the forward flight condition. The 
distribution in the multiobjective formulation case shows the trade-offs associated with these two 
conflicting requirements. From root to mid span, the distribution is much closer to the airplane 
mode twist, however, at outboard locations the optimal twist is nearly identical to the hover only 
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case. The reason the multiobjective optimal twist distribution more closely follows the cruise only 
case is due to the fact that the rotor planform is sized for hover and therefore in order to 
simultaneously achieve high propulsive efficiency it is necessary to twist the blade according to the 
requirements in cruise. It is noteworthy that the optimal twist distribution in case of the 
multiobjective formulation is almost identical to the reference rotor. 



Figure 7 Chord distributions 


Figure 9 shows that the zero lift angle of (oc z i) attack is drastically reduced (made more 
negative) throughout the blade span. This effectively increases the camber of the airfoil which 
improves the lift-to-drag ratio and improves performance. In cruise, there is an increase in a z i 
from the reference value at the blade tip. The reason for the resulting decrease in camber at the tip 
can be explained as the optimizer’s effort in avoiding a reduction of the drag divergence Mach 
number (Mdd) to values below operational Mach number at this location. Once again, the 
multiobjective formulation more closely follows the hover distribution except at radial locations 
which require high Mach drag divergence number (e.g. near the tip). The results clearly 
demonstrate the trade-offs associated between the two flight conditions. It must be noted that the 
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zero lift angle of attack distributions presented in this paper are achieved primarily through blade 
camber and a more comprehensive study would have to include considerations of not only 
maximum camber requirements but must also include constraints on the moment coefficients as 
well. 



Nondimensional radius, y/R 


Figure 8 Twist distributions 


Figure 10 shows the thickness to chord distributions (t/c) where it is seen that the blade 
thickness is significantly decreased in all cases, except at the blade tip. These decreases are an 
attempt to improve the rotor performance by decreasing the profile drag through reductions in the 
airfoil thickness. It is of interest to note that the multiobjective case is heavily biased towards the 
hover solution. This is explained as follows, since no weight factors are used in the K-S objective 
function formulation and all of the objective functions and constraints are normalized to unitary 
values, the optimization procedure is driven primarily by the dominance of the hover thrust 
requirement. 
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Reference 



Nondimensional radius, y/R 


Figure 9 Zero lift angle of attack distributions 



Nondimensional radius, y/R 


Figure 10 Thickness to chord ratio distributions 
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The final optimization case is performed for the simultaneous maximization of the hover figure 
of merit and propulsive efficiency of the XV- 15 rotor at a forward speed of 400 knots. Results are 
compared with a reference rotor with the same planform as used in the 300 knots case and also the 
optimum rotor obtained from the multiobjective function formulation case at 300 knots. In Fig. 1 1 
it is shown that the increase of 5.5 percent in the hover figure of merit, in case of the 400 knots 
cruise condition, is not as large as in case of the 300 knots cruise condition (7.7 percent). It is 
further shown that although the blade was originally designed for 300 knots, the same level of 
efficiency can be achieved at 400 knots after optimization (0.5 percent increase). The reason this 
increase in the figure of merit, from reference to optimum, is not as large in the case of the 400 
knots cruise condition is due to the more stringent requirements of the cruise condition. 

| | Reference (300 Knots) 
p/j Optimum (300 Knots) 

I H Optimum (400 Knots) 

3.7% 



Figure 1 1 Comparison of optimum results 

The chord distributions for the reference and the optimum rotors for the cruise speed of 400 
knots are presented in Fig. 12. The multiobjective case with a cruise speed of 300 knots is 
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presented for comparison. The figure shows significant changes in chord from 300 to 400 knots 
designs. In case of the 400 knots forward flight condition, the chord is increased from midspan to 
outboard. The drastic reduction in chord at the root is due to an inability of the optimizer to find a 
satisfactory, that is, unstalled solution consistent with the cruise requirements. 

The effective chordwise angle of attack (a) distribution in hover is shown in Fig. 13 where 

a(y) = e(y)+e 7 5-«i(y)- a zi(y) ( 10 > 

and 075 is the collective pitch. From the figure it is seen that the effective angle of attack at the root 
in case of the reference blade and the 400 knot multiobjective optimization case is much higher than 
the effective angle of attack of either the hover only optimization case or the 300 knot multi- 
objective optimization case. 

— — Reference 



Nondimensional radius, y/R 
Figure 12 Chord distributions 


The sectional lift-to-drag coefficient ratio (q/cd) is presented in Fig. 14 and indicates that this 
large angle of attack at the root causes blade stall. The optimizer attempts to overcome this problem 
by reducing the root chord in an effort to decrease the stall drag penalty. The decreased values of 
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the chord, from reference to optimum, in the 400 knot multiobjective case can be explained by 
examining Figs. 15 and 16. It is seen in Fig. 15 that in case of the reference blade at 400 knots, 
the effective chordwise angle of attack (a) in cruise is very small over the midspan of the blade and 
therefore this portion of the blade provides very little lift (Fig. 16). By reducing the chord in the 
optimum blade over the midspan region, the angle of attack is redistributed more uniformly 
throughout the blade. As a result the lift is more evenly distributed and the performance is 
improved. 



Figure 13 Effective chordwise angle of attack in hover 


The twist distributions are shown in Fig. 17. From the figure it is seen that the optimum twist 
distribution for the 400 knots case has slightly less total twist than the reference rotor. The 
distribution is also more linear, compared to the reference. This reduced twist and linearity result 
from the higher axial advance ratio. To generate more thrust from the midspan sections, the angle 
of attack is increased slightly in this region. This has the effect of making the overall twist 
distributions more linear. Again, the similarities in the distributions between the reference and the 
300 knots optimum rotor must be noted. 
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Figure 16 Blade thrust distribution in cruise (400 knots) 


Reference 


O 300 Knots 



Nondimensional radius, y/R 

Figure 17 Twist distributions 
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Figure 18 shows the zero lift angle of attack (a z j) distributions. It is seen that the optimizer 
dramatically reduces a z i (increasing camber) over bulk of the blade span returning to the reference 
value at the tip. In comparing the two cruise cases, over the outboard section, there is significant 
reduction in camber in the 400 knots design compared to the 300 knots design. This is because at 
400 knots the rotor is operating at a relatively higher effective Mach number than it is at 300 knots. 
Therefore, the reduction in camber is necessary to ensure that the blade sections of the 400 knots 
rotor operate below Mdd to avoid large drag penalties. 



1 1 1 1 1 1 1 1 I I 

0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 

Nondimensional radius, y/R 

Figure 18 Zero lift angle of attack distributions 


The thickness to chord ratio (t/c) distributions for the reference and the optimum rotors are 
shown in Fig. 19. The distributions are similar in both the 300 and the 400 knots cases at inboard 
and midspan sections. In both cases, significant reductions occur throughout the blade span. 
Such reductions in the t/c ratio are the result of an effort to increase Mdd which is accomplished 
primarily through thinner airfoils. The root value of t/c is slightly higher in case of the 400 knots 
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rotor, compared to the rotor with a cruise speed of 300 knots. This anomaly can only be explained 
by noting that since the lift is so small at this location for both hover and cruise (Figs. 14 and 16), 
sensitivities to changes at the root are very small and therefore their effect on rotor performance in 
either flight condition is minimal. The tip values of the thickness to chord ratio are slightly 
decreased in case of the rotor operating at 400 knots forward speed, compared to the 300 knots 
case, due to the fact that the drag divergence Mach number must be significantly increased at these 
locations to ensure that the they are above the elemental Mach number. 



Figure 19 Thickness distributions 

Off-design studies for the hover figure of merit and the high speed cruise propulsive efficiency 
are shown in Figs. 20 and 21. Figure 20 illustrates that the hover figure of merit remains 
reasonably high over small increases in the coefficient of thrust (Ct). The performance is better in 
the 400 knots multiobjective design case. It is interesting to note from Fig. 21 that the 
multiobjective designs at both 300 and 400 knots maintain reasonable values of ri c at high off- 
design values of Cy. The efficiency drops down drastically at higher Cy values in the 300 knots 
cruise-only optimization case. This can be explained by the fact that in this case, the propulsive 
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efficiency is increased during optimization by trimming to values very near blade stall in an effort 
to increase the lift-to-drag ratio of the airfoil. As a result, as the thrust is increased, a greater 
collective pitch is required which causes the blade to stall. This is caused by the unrealistic chord 
distribution near the root in this case. 

To summarize, the procedure developed shows the adequacy of the optimization technique in 
providing important design trends and trade-offs associated with improving aerodynamic efficiency 
in cruise and in hover sequentially as well as simultaneously. The procedure will serve as a first 
level in a multilevel optimization procedure. Therefore, the numerical results presented in this 
paper, purely obtained using performance criteria, are expected to change as other disciplines are 
integrated. For example, the high camber which results in high LTD ratios can also lead to high 
pitching moments. This problem can be addressed by combining structural and dynamic design 
criteria in the optimization procedure. 



Figure 20 Hover figure of merit off-design study 
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Figure 21 High speed cruise propulsive efficiency off-design study 


3.7 Concluding Remarks 

A performance based optimization procedure is developed to investigate the selection of blade 
aerodynamic characteristics for improving propulsive efficiency in high speed cruise as well as the 
figure of merit in hover. The classical blade element momentum approach is used for the 
aerodynamic analysis. The code includes an empirical correction to two dimensional stall behavior 
in order to represent the high lift capability of rotors and propellers. The computational ease 
offered by the code allows for several design optimization problems to be studied and the 
subsequent trade-offs investigated. Several different flight conditions and optimization problems 
are considered. The following important observations are made. 

1 . Maximum improvements in aerodynamic performance were obtained from the single objective 
function optimization procedures. The optimization algorithm converged to the well known 
classical optimum solutions for hover and propellers. The multiobjective function case 
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demonstrated the trade-offs associated with the conflicting requirements of the two flight 
conditions. 

2. Improvements were made in the hover figure of merit and cruise propulsive efficiency in all 
cases and at both 300 and 400 knots. The improvements at 400 knots were not as significant 
due to the more demanding cruise requirements at this speed. The aerodynamic performance in 
this case was similar to the performance obtained in the 300 knot case, which was achieved 
through the proper tailoring of the blade spanwise properties. 

3. The changes in the thickness and camber distributions of the optimum rotors represented an 
effort to increase the drag divergence Mach number without degradation of the blade lifting 
characteristics. The rotor operating at a forward speed of 400 knots required less camber than 
the 300 knots rotor in order to achieve the higher drag divergence Mach numbers necessary to 
remain above elemental Mach numbers. 
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4. Coupled Aerodynamic and Structures Optimization 


For this objective, the aerodynamic and structural design criteria in both high speed cruise and 
in hover are addressed by developing a multilevel decomposition based optimization procedure. At 
the upper level, the aerodynamic performance of prop-rotors is optimized for both high speed 
cruise and hover using planform variables. Constraints are imposed on the rotor thrust in hover 
and in cruise. At the lower level, the rotor is optimized for improved structural performance using 
composite ply orientations as design variables. Since the lower level optimization involves discrete 
design variables, an optimization procedure, based on simulated annealing technique [14], is 
developed to address this complex problem. The results of the optimization procedure are 
compared with the same reference rotor used in the first objective (Chapter 3). 

For the aerodynamics analysis, the same algorithm used in the aerodynamic optimization 
objective (Chapter 3) is used. The procedure offers a significant reduction in computational effort 
from more comprehensive procedures which were previously used by Chattopadhyay and 
McCarthy [4,5,15-17] as demonstrated in the first objective of this report. The use of this analysis 
within an optimization procedure also provides realistic design trends as demonstrated by 
McCarthy et al. [18]. The structural analysis is performed using a composite box beam model 
which includes blade pre-twist, taper and sweep to represent the principal load carrying member in 
the blade. In this research, a quasi one-dimensional composite beam theory, based on the model 
developed by Smith and Chopra [19], is extended to include beam pre-twist, taper and spanwise 
sweep. The procedure offers significant computational savings from more refined finite element 
models while maintaining sufficient accuracy. 

4. 1 Multilevel Optimization Problem 

Multidisciplinary design optimization of rotary wing aircraft can be a computationally intensive 
task if all of the disciplines involved are coupled with the optimizer and the optimization is 
performed “all-at-once” in a single level. Decomposition techniques, through which such complex 
optimization problems can be reduced into a number of sub problems, can be very effective in 
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addressing this issue. Very recently, such techniques have been applied to rotary wing 
optimization problems [20,21], In this paper, the multilevel decomposition procedure is used to 
address the complex multidisciplinary issues associated with high speed prop-rotor design. The 
optimization procedure is decomposed into two levels. At the upper level the aerodynamic 
performance is improved and at the lower level, the objective is to improve the structural 
performance of the blade using composite tailoring. The upper and lower levels are coupled 
through the use of optimal sensitivity parameters [21], which are essential in maintaining proper 
coupling between the levels. Following is a description of the decomposition and the optimization 

problem formulation. 

4.1.1 Upper Level 

At this level, the axial efficiency in high speed cruise (Tfc) and the hover figure of merit (FM) 
are maximized simultaneously. Constraints are imposed on the physical dimensions of the blade to 
ensure that the load carrying member of the rotor is maintained within the dimensions of the airfoil. 
The blade is discretized and design variables include the values of the chord (c), twist (9), 
thickness to chord ratio (t/c), zero lift angle of attack «X z i) at each node. The sweep distribution, 
however, is not discretized in order to ensure continuity of the elastic axis. The lifting line is 
assumed to be a quadratic function of the following form. 

Xa/c = dix + d 2 (4y 2 - 3y) (10) 

where Xa/ C is the position of the aerodynamic center measured from the leading edge of the chord. 
The aerodynamic center is assumed to coincide with the elastic axis of the load carrying structural 
member used in the problem formulation. The coefficients which determine the position of the 
aerodynamic center, d, and d 2 , are used as design variables. Note that the functions associated 
with these coefficients (f,(y) = y and f 2 (y) = 4y* - 3y) are orthogonal to each other which is a 
favorable characteristic for optimization. The sweep distribution can then be formulated using this 

expression as follows. 
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( 11 ) 


A(y) = tarr^dx/dy) 

= tarHfdi +d2(8y - 3)} 

The optimization problem is stated as follows. 

Maximize FM(<I>), T| c (d>) 


subject to T h = (Th) re f 
Tc = (T c ) re f 
Kthor — Imax 

where O = [c(y), 0(y), A(y), a z i(y), t/c(y)] is the design variable vector. The quantity th or is the 
thickness of the horizontal wall in the box beam, t ma x is the maximum thickness of the airfoil and jc 
is a scaling factor used to ensure that the box beam is maintained within the airfoil cross section. 
Since the optimization problem involves more than one design objective, the problem is formulated 
using the Kreisselmeier-Steinhauser (K-S) function approach (Chapter 3). 


4.1.2 Lower Level 

The objective at this level is to minimize the tip displacements in hover and in cruise. The most 
critical displacements in hover are the vertical displacement (wh) and the elastic twist (<|>h). In 
cruise, the elastic twist (4> c ) and the inplane displacement (v c ) are important. Therefore, these four 
displacements are selected as the individual objective functions to be minimized. Ply orientations 
are used as design variables. However, to avoid impractical orientations, the ply angles are chosen 
from a set of standard values [0°, ±15°, ±30°, ... , 90°]. The Tsai-Wu failure criterion [22] is used, 
which assumes that to avoid material failure the following equation representing a failure surface in 
stress-space must be satisfied. 



where Gi and <?2 represent normal stresses along the material axes and X \2 represents the shear 
stress (see Fig. 22). The subscripts (T), (C) and (s) represent the ultimate stress in tension, 
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compression and shear, respectively. This reduces the total number of constraints as individual 
constraints on the stress (a i, 02 and 1 12 ) at each ply are avoided. Each of the composite plates 
used in the box beam modeling are assumed to be symmetric about the midplane of the plate and 
the beam itself is assumed to be symmetric about its local axes, r\ and C, (Fig. 23). Therefore, the 
above failure criterion is imposed, on each lamina, at each of the four comers of the box beam to 
prevent failure due to stresses. 



Figure 22 Composite lamina material axes 



Figure 23 Blade cross section 


4.2 Analysis 


4.2.1 Aerodynamic Analysis 

The aerodynamic formulation is based on the model used in the aerodynamic optimization 
objective. Further details of this analytical technique are found in Chapter 3 and Ref. 18. The 
aerodynamic analysis is coupled with a structural analysis which is described next. Therefore, the 
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results of the complete analysis represent trimmed static solutions which include the effects of 
elastic deformations. 



Figure 24 Composite box beam 


4.2.2 Structural Analysis 

The load carrying member of the rotor is modeled as a single cell composite box beam 
(Figs. 23 and 24). In addition to the structural member inside the airfoil, the weights of the 
honeycomb structure and the blade skin are also included in the total weight calculation. The blade 
is discretized using Finite elements with 19 total degrees of freedom, U e , and unequal element 
sizes. Using the finite element model it is possible to incorporate blade pre-twist and sweep 
distributions into the problem formulation. The nodal degrees of freedom are described as 
follows. 

Ue T = [ui,U2,U3,U4,Vb 1 ,Vb 1 ,Vb2,Vb2» w bi>Wbi,Wb2»Wb2,<l)l,<l)2.<l>3.Vsi,Vs2,W sl ,W S 2) (13) 

where u is the axial displacement, v and w are the inplane horizontal and vertical displacements and 
<j> is the elastic twist. First partial derivatives with respect to the spanwise axis (x) are denoted ( ' ). 
The formulation assumes that the inplane displacements can be decoupled into a term 
corresponding to pure bending and a term corresponding to shear as shown below. 
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V = V b + v s 


(14) 

(15) 


W = w b + w s 


where the subscript (b) refers to the displacement due to beam bending and the subscript (s) refers 
to the displacements due to shear. Identical node locations are used in specifying both 
aerodynamic and structural parameters. 

The outer dimensions of the box beam (Fig. 24) are constant percentages of the chord and 
thickness. Each composite plate used to model the composite box beam is assumed to be 
symmetric about the midplane of the plate. Each plate is made up of 24 laminated orthotropic 
composite plies. Further, the box beam is assumed to have double symmetry about the local 
coordinate axes (r| and Q. This ensures that the two vertical walls are identical to each other and 
also that the two horizontal walls are identical. The vertical and horizontal walls, however, are 
assumed to be independent of each other. The beam cross section is described by stretching, 
bending, twisting, shearing and torsion related warping. The stress-strain (a - e) relationships in 
the vertical and horizontal walls are written as follows. 
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(17) 


where the subscripts (v) and (h) indicate vertical and horizontal walls, respectively. The off-axis 
stiffness matrix for each lamina is denoted Qjj (i, j = 1,2 & 6). The strains in the vertical walls are 


expressed as follows. 


e xx - u ' - T|(v" * YxV) ' 

(18) 

( 5/\. w ) 


= r\- V + y ft 

V ^ J 

(19) 
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aou' + bir|(v" - y^n') - c 2 £(w" - y&') - (d G + dm + d 2 Q<j>' 
+ (f 0 + Wy xh + (go + gi)Tx"; 


( 20 ) 


£ cc = 


where u is the axial displacement, v and w are the horizontal and vertical inplane displacements, 
respectively, <() is the twist angle and y^ and yj£ represent the inplane shear stresses (Fig. 24). 
The superscript ( " ) represents second partial derivatives with respect to x. Similar expressions 
are obtained for the horizontal walls. The out-of-plane warping is denoted X w and the coefficients 
ao - g] are determined such that the net inplane forces and moments in the horizontal and vertical 
walls are zero. Using these conditions, the equations of equilibrium are written as follows. 
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where M z and M y are the lagging and flapping moments, respectively and T is the torsional 
moment. The axial force is denoted F, the inplane horizontal and vertical shear forces are denoted 
Q y and Q z , respectively and ky (i, j = 1,2, ... , 6) represent the stiffness matrix elements. Further 
details of the formulation can be found in Ref. [19]. 

The elemental equilibrium equations are written in vector form as follows. 


Fe — KgUg (23) 

where F e = [F Q y Q z T m x m z ] T is the elemental force vector, K e = [kjj] is the elemental 
stiffness matrix and u e = [u' y^ Jyfc <}>' (v" - y^') (w" - y^')] T represents the elemental degrees 

of freedom in the local coordinate system. Using these principles a finite element approach is 
developed using the weak formulation as follows. 
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( 24 ) 


R R 

J H G F G dx = J H G K G H G U G dx 

O O 

where Fg and Kg are the global force vector and stiffness matrix, respectively, Ug is the global 
vector of discrete nodal degrees of freedom and Hg is the global representation of the interpolation 
matrix. The blade displacements are written in the discretized elemental form as shown below. 

u(s) = H(s)U e (25) 

where H(s) represents the complete set of interpolation vectors and U e is the elemental nodal 
displacement vector (Eqn. 13). 

4.3 Optimization Implementation 

The multilevel optimization technique used in this research decomposes the design problem into 
a level where only continuous design variables are used (upper level) and a second level where 
discrete design variables are used (lower level). Therefore, two separate optimization strategies are 
necessary. Since traditional optimization techniques minimize only a single objective function, it is 
necessary to transform the original objective functions, at each level, using multiobjective function 
techniques. The technique used for each level is the Kreisselmeier-Steinhauser (K-S) function 
approach. A nonlinear programming procedure (NLP) based on the Davidon-Fletcher-Powell 
method (DFP) is used as the optimization algorithm at this level [23]. For the lower level, an 
optimization procedure has been developed based on the method of simulated annealing [10] which 
is described next. 

If ply angles are used as continuous design variables and are later rounded off to the nearest 
practical value (e.g. 47.3° being rounded off to 45°), the result can lead to sub optimal designs 
[24]. Allowing the ply angles to vary within a set of prescribed values during optimization is a 
more efficient approach. Therefore, in the structural optimization problem, the plies are selected 
from within a set of pre-selected orientations of integer multiples of +15° such as 0°, ±15°, ±30°, 
..., 90°. This reduces the problem into a discrete optimization problem and conventional gradient 
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based optimization procedures are no longer applicable. Using the procedure developed in the 
current work, the optimization problem is transformed into a completely discrete problem. The 
new discrete problem can now be optimized using a simulated annealing technique which has been 
shown to be effective in a variety of different engineering applications [25-27], The procedure 
developed is outlined below. 

The multiobjective optimization problem is stated as follows: 

Minimize F K (aj) K = 1, 2, NF 

subject to gM(aj) M = 1, 2, NC 

i = 1,2, NDV 

where NF is the number of objective functions, NC is the number of constraints, NDV is the 
number of design variables and a is the vector of design variables. The i th design variable can 
assume any value from the design variables vector of pre-selected values, dj q such as 0°, ±15°, 
±30°, ... , 90° for ply angles. Therefore, the objective functions and constraints are completely 
represented by discrete design variables. The multiobjective optimization problem is transformed 
into a single unconstrained composite function to be minimized using the K-S function approach as 
previously described. An optimal solution cannot be guaranteed for a discrete optimization 
problem without evaluating every possible combination of discrete variables which is 
computationally impractical. Near optimal solutions can be obtained, however, with significant 
improvements in all objective functions with a reasonable amount of computational effort. 

The simulated annealing algorithm is described briefly as follows. 

START 

Current design is F 

Perturb current design F new 

V Fnew < F then 
F = F new 

Else if p a cc £ p then 
F = F n ew 

End if 

Go to START 
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where F is the objective function to be minimized and p is a random number such that 0 < p < 1. 
The acceptance probability p acc of retaining a worse design is computed as follows. 

Pace = ex p(4f ) (26) 

where AF is the change in the objective function and T k is the “temperature” which is reduced 
during successive iterations, to ensure smooth convergence, according to: 

T k = (r c ) k T 0 (27) 

where T 0 is the initial temperature and r c is the cooling rate which determines the temperature at the 

k th iteration. This reduces the probability of accepting a worse design. Occasionally accepting a 
worse design under the given probability allows the algorithm to climb out of possible local 
minima. The above loop is repeated a prescribed number of times for each cycle in the multilevel 
optimization procedure. 

4.4 Results 

The reference rotor used is the same representation of the XV- 15 prop-rotor which was used in 
the aerodynamic optimization objective. Identical flight conditions are used which are summarized 
in Table 1. The blade is discretized into 10 segments (11 node points). The composite material 

used in the structural analysis is carbon-PEEK AS4/APC2 [28] (Table 2). 

At the upper level the design variables for the chord (c), twist (0), zero angle of attack (a z i) 
and thickness to chord ratio (t/c) are all based on discrete nodal values whereas the sweep 
distribution (A) is assumed to be based on a quadratic lifting line. This yields a total of 46 design 
variables. The scaling factor (k), used in the upper level constraints to ensure that the box beam is 
maintained in the airfoil, is assumed to vary as follows. 


'4.0 

y <0.8 


3.0 

0.8 <y <0.95 

(28) 

2.5 
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0.95 < y < 1.0 
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where y is the nondimensional radial location. The reason for this variation is to ensure the 
validity of thin wall theory at inboard sections of the blade, since the majority of the load is carried 
in this section, without being too restrictive at the tip where thinner airfoils can improve the 
aerodynamic performance. The tolerance on the lower level objective functions, used as optimal 
sensitivity parameters at this level, is initially 5 percent. This value is reduced to 0.5 percent after 6 
complete multilevel cycles in order to ensure that there is very little degradation of lower level 
objective functions at the upper level. A nonlinear programming technique based on the DFP 
algorithm, developed at Arizona State University, is used for optimization at this level. 

Table 2 Composite properties 


Carbon-PEEK AS4/APC2 


Ei 

19.43 

(106) 

p.s.i. 

e 2 

1.29 

(106) 

p.s.i. 

Gj2 

0.74 

(106) 

p.s.L 

Vl2 

0.28 



(Jit 

309. 

(103) 

p.s.i. 

<*1C 

160. 

(103) 

p.s.i. 

02j 

11.6 

(103) 

p.s.i. 

Ct2C 

29.0 

(10 3 ) 

p.s.i. 

^12 s 

23.2 

(lO 3 ) 

p.s.i. 


At the lower level, the discrete values of the ply orientations are used as design variables. 
Since a symmetric and balanced lay-up is used in both the vertical and horizontal walls, this leads 
to 12 independent design variables which can assume any of the 7 pre-selected values of ply angle 
orientations. The design space associated with this problem therefore consists of nearly 14 billion 
(7 12 ) possible combinations. The actual number of loops used at each multilevel cycle is 1000 
iterations which represents only a very small portion of the total design space. The value of the 
K-S function multiplier (8) used at this level is initially set to 5.0, but is subsequently reduced to 
1 .0 after 4 multilevel iterations. This is done to increase the sensitivity of the K-S function to 
changes in all of the objective functions and constraints. The cooling rate (r c ) used in the simulated 
annealing algorithm is 0.995 and the initial temperature (To) is set to 1.0. 
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Convergence of the multilevel optimization procedure is achieved in 7 cycles, where a cycle 
represents an individually converged optimization including both upper and lower levels. The 
initial design used in this optimization procedure corresponds to the optimal results obtained from a 
purely aerodynamic optimization using the same aerodynamic analysis (Chapter 3) [18]. The 
results from the multilevel based optimization are compared to the design obtained from the purely 
aerodynamic optimization and also to a reference rotor which represents the original XV- 15 prop- 

rotor. The results are presented in Table 3 and Figs. 25-33. 

The upper level objective functions are presented in Fig. 25. It is shown in the figure that the 
hover figure of merit (FM), which is increased by 6.6 percent in the purely aerodynamic 
procedure, is further increased in the multilevel optimization (7.5 percent) compared to the 
reference rotor. The propulsive efficiency in high speed cruise Cn c ) is increased by 3.6 percent 
from the reference rotor in the aerodynamic formulation and in the multilevel procedure the increase 
is slightly less (3.2 percent). These trends can be explained by examining the rotor planform. 

□ Reference 
^ Aerodynamic 
^ Multilevel 



FM ti 

Figure 25 Summary of upper level objective functions 
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The optimum and reference chord distributions (c) are shown in Fig. 26 where it is observed 
that the planforms obtained from both the aerodynamic and the multilevel optimization procedures 
are similar in shape and they differ sharply from the reference chord distribution. The optimum 
distributions closely resemble the well known ideal hover planform. Such a distribution is the 
result of the large solidity requirement necessary to achieve the required thrust in hover. The 
reduction in the outboard section, relative to the reference rotor, is attributable to the lack of any 
maneuver margin requirement in the optimization problem formulation. 
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Figure 27 shows the thickness to chord ratio (t/c) distributions for the three rotor 
configurations. From the figure, it is seen that there are large reductions from the reference rotor at 
the inboard sections in both of the optimum rotors. In case of the purely aerodynamic formulation, 
the optimum thickness is reduced from the reference blade throughout the blade span except at the 
90 percent span location. This decrease represents an attempt to improve the rotor performance by 
reducing the profile drag. In case of the multilevel optimization procedure, similar trends are 
noted, only the reductions from the reference rotor are not as large. In fact, at the outboard section 
the airfoil thickness is actually increased from the reference values. The reason for these increases 
at the tip are to satisfy the constraint which ensures that the box beam can be contained within the 
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airfoil. In case of the aerodynamic procedure, these constraints are violated at the tip. Since the 
multilevel procedure results in a thicker blade, compared to that from the aerodynamic 
optimization, the maximum coefficient of lift (CL max ) is increased in this design. This increase in 
C Lmax subsequently improves the hover figure of merit which explains the slight increase in FM 
from aerodynamic formulation. However, associated with the thicker airfoils are lower drag 
divergence Mach numbers (Mdd) which cause the high speed cruise performance to be slightly 
diminished, from the purely aerodynamic formulation. 

Reference 



Figure 27 Thickness to chord ratio distributions 

The lifting line (Xa/ C ) and sweep (A) distributions for the multilevel optimization and the 
reference rotor are presented in Figs. 28 and 29. Note that in the purely aerodynamic formulation 
no sweep modeling is included. As is shown in Fig. 28, the lifting line is swept slightly forward, 
with the maximum offset at the tip. This results in the near linear distribution of the sweep. This 
sweep distribution has the effect of increasing the torsional moment about the hub in an effort to 

reduce the elastic twist, which in this case is negative (nose down). It must be noted that to 
investigate the effect of forward sweep, which can have an adverse affect on aeroelastic stability, a 
more comprehensive formulation including dynamic loads would have to be used. 
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The twist (9) and the zero lift angle of attack (a z j) distributions are shown in Figs. 30 and 31. 
From Fig. 30 it is seen that the twist distribution for the aerodynamic and the multilevel 
optimization results are nearly identical varying only slightly from the reference distribution. The 
variations are most notable at the midspan locations where the optimum twist, through this region, 
is more linear than the reference rotor. This is due to the high speed cruise requirements where 
more linear twist distributions are expected due to the high inflow velocity in forward flight. 
Figure 3 1 shows that the a z | distributions for both the aerodynamic and the multilevel optimization 
procedures are nearly identical. However, in this case, both of these distributions vary drastically 
from the reference distribution. The figure shows that the optimum distributions are decreased 
significantly from the reference values at the inboard and midspan locations. This effectively 
increases the camber of the airfoil which improves the airfoil lift-to-drag ratio thereby improving 
performance. The zero lift angle of attack is then reduced (less camber) at the outboard section and 
is increased at the blade tip. This can be explained as the optimizer’s attempt to avoid the reduction 
of Mdd to values below the operational Mach number at this location. 
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Figure 28 Lifting line distributions 
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Figure 29 Blade sweep distributions 



Figure 30 Blade twist distributions 
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Reference 



Figure 3 1 Zero lift angle of attack distributions 

The nondimensional values of the lower level objective functions are presented in Fig. 32 
where large deviations from the reference rotor are observed in the purely aerodynamic and 
multilevel optimization procedures. It should be noted here that the purely aerodynamic results 
were obtained in Ref. 18 (Chapter 3) without any structural considerations. The structural 
performance results pertaining to that case, which are presented here, are computed using the 
optimized rotor planform obtained from that study. In the aerodynamic formulation it is seen that 
the elastic twist in hover (<J>h) and in cruise (<j> c ) and the vertical displacement in hover (w^) are 
significantly increased by 53, 61 and 6.5 percent, respectively, from the reference rotor. 
However, the horizontal displacement in cruise is actually decreased by 15 percent. In case of the 
multilevel optimization in which these displacements are included as objective functions, it is seen 
that there are significant reductions in all four displacements, from reference to optimum, after 
optimization. In this case, the elastic twist in hover and in cruise are reduced by 9.9 and 14 
percent, respectively and the vertical displacement in hover and the horizontal displacement in 
cruise are reduce by 16 and 3.9 percent, respectively. The magnitude of all the displacements are 
also reduced throughout the blade span from the reference values. An example of this is seen in 
Fig. 33 which displays the elastic twist deformations. This shows that in case of the multilevel 
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formulation, even though the aerodynamic performance is comparable to the purely aerodynamic 
design (Fig. 25), the structural performance is greatly improved. These reductions are attributable 
to increased thickness of the airfoil, which in turn increases the height of the box beam, and also to 
the improved stacking sequence. 

□ Reference 
^ Aerodynamic 
^ Multilevel 





<t>h Wh 4>C v c 

Figure 32 Summary of lower level objective functions 

The composite laminate stacking sequences are presented in Table 3. Note that in the reference 
blade, the horizontal and vertical walls are assumed to have the same stacking sequence. Further, 
since all of the laminates are considered to be symmetric about their midplane, only 1 2 of the 24 
total plies are presented in the table for each wall. The individual ply thickness used in this study is 
0.001 inches which results in a total wall thickness of 0.24 in. For the sake of comparison, the 
structural response in case of the purely aerodynamic optimization formulation is obtained using 
the reference stacking sequence. Note that in the optimal configuration, there is a reduced number 
of +45° plies and an increased number of ±30° plies. This can be explained as an attempt to 
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compromise between the conflicting requirements of reduced elastic twist and reduced transverse 
displacements. 


Reference (hover) 

■Q — Optimum (hover) 


Reference (cruise) 



Figure 33 Elastic twist distributions 


Table 3 Ply orientation angles 


Reference Optimum 

Horizontal Wall Vertical Wall 


Outer ply CP 15° CP 


— 

CP 

-15° 

CP 


CP 

CP 

15° 


CP 

CP 

-15° 


15° 

15° 

30P 


-15° 

-15° 

-30° 


15° 

15° 

15° 


-15° 

-15° 

-15° 


45° 

30P 

45° 


■ 45 ° 

-30° 

^45° 


45° 

30P 

15° 

— 

midplane -45° 

-30° 

-15° 
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4.5 Concluding Remarks 

A multilevel optimization procedure based on aerodynamic and structural performance is 
developed for an investigation of the coupled aerodynamic and structural design problem of high 
speed prop-rotors. The classical blade element momentum approach is used for the aerodynamic 
analysis. The structural analysis is performed using a quasi one-dimensional finite element 
approach based on a composite box beam model. At the upper level, a nonlinear programming 
technique based on the Davidon-Fletcher-Powell method is used as the optimization algorithm. At 
the lower level, a discrete optimization based on the simulated annealing technique is used. The 
optimization problems are formulated using the Kreisselmeier-Steinhauser function approach. The 
optimum results obtained in this study are compared with a reference rotor and with a design 
obtained using a purely aerodynamic optimization procedure. Convergence is achieved after 147 
total aerodynamic optimization cycles at the upper level and a total of 7000 simulated annealing 
iterations at the lower level. The following important observations are observed. 

1 . The multilevel optimization procedure significantly improves the aerodynamic and structural 
response of a high speed prop-rotor blade. Further, the aerodynamic performance is 
comparable to a purely aerodynamic design demonstrating the success of the decomposition 
based procedure. 

2 . The simulated annealing algorithm successfully minimizes the tip displacements by altering the 
composite plate stacking sequences in the horizontal and the vertical walls. The optimum 
composite stacking sequence represents a compromise between reduced elastic twist and 
reduced transverse bending. This is manifested through the selection of ±30 degree plies in 
both the horizontal and the vertical walls. 

3 . The airfoil thickness to chord ratio is increased from the purely aerodynamic formulation due to 
structural constraints on the blade. The thicker airfoil increases the hover performance with 
only a slight degradation in the cruise efficiency from the optimum aerodynamic results. 
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4. The optimum blade is swept slightly forward to improve the elastic twist. The sweep increases 
the torsional moment which has the effect of reducing the nose down elastic twist. 

5. The chord, twist and zero lift angle of attack distributions differ significantly from the reference 
blade but are very similar to the optimum aerodynamic configuration. 

6. The elastic deformations at the blade tip as well as throughout the blade span are reduced from 
reference values. This is achieved through a combination of an improved ply stacking 
sequence and slightly thicker airfoils. 


46 



5. Multipoint Optimization 

The optimization procedure of Chapter 4 is now extended to include the take-off performance 
in the formulation. The aerodynamic and structural design criteria in high speed cruise, in hover 
and in take-off are addressed using the multilevel decomposition-based optimization procedure 
developed in Chapter 4. At the upper level, the aerodynamic performance of prop-rotors is 
optimized for each flight condition using planform variables. Constraints are imposed on the rotor 
thrust in all three flight conditions. A nonlinear programming technique based on the Broyden- 
Fletcher-Goldfarb-Shanno (BFGS) [23] algorithm is used for the optimization. At the lower level, 
the rotor is optimized for improved structural performance using composite ply stacking sequence 
as design variables. Since only discrete design variables are used at the lower level, an 
optimization procedure, based on simulated annealing algorithm (Chapter 4), is developed to 
address this complex problem. 

5.1 Analysis 

The aerodynamic and structural analyses are performed in this objective with the same 
formulation used in the Chapters 3 and 4. Details of these formulations can be found in those 
chapters. 

5.2 Optimization Implementation 

The primary objectives of this study are to improve aerodynamic and structural performance of 
the tilt-rotor aircraft at three design points, high speed cruise, hover and take-off. Since the 
problem is complex and is associated with several objectives functions, constraints and design 
variables, a multilevel decomposition technique is used to decompose the problem into two levels. 
In this study, the decomposition technique developed in Chapter 4 (Ref. 29) is used to formulate 
an optimization problem involving a total of three flight conditions. The optimization procedure is 
decomposed into two levels. The aerodynamic performance is improved at the upper level and the 
structural criteria are addressed at the lower level. Following is a description of the two levels. 


47 



5.2.1 Upper Level 

The axial efficiency in high speed cruise (T| c ) and the figure of merit in both hover and in take- 
off (FMh and FM t , respectively) are maximized simultaneously in this level using aerodynamic 
design variables. Constraints are imposed on the rotor thrust at each of these flight conditions. 
Geometric constraints are also imposed on the physical dimensions of the blade to ensure that the 
load carrying member of the rotor is maintained within the dimensions of the airfoil. The blade is 
discretized and design variables include the values of chord (c), twist (0), thickness to chord ratio 
(t/c) and zero lift angle of attack (a z i) at each node. To ensure monotony of the sweep, the 
following quadratic variation is used to represent the lifting line. 

Xa/c = 5 d iy 2 ( ' 29 ' ) 

where ya/ c is the position of the aerodynamic center, which in this formulation coincides with the 
shear center. The coefficient di is used as a design variable to determine this position. Note that 
this sweep formulation is slightly different from the formulation used in Chapter 4. The sweep 
distribution can then be formulated using this expression as follows. 

A(y) = tan _1 (diy) (^0) 

This also ensures continuity of the elastic axis. The optimization problem is stated as follows. 
Maximize 

tlc(<D), FM h (<D), FM t (0) 

subject to 

T c = (Tc)ref 

Th = (Th)ref 

T, = (Tt) re f 

Kthor — tmax 

where the subscript ‘t’ refers to the take-off condition and the other subscripts are as described 
before. 
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5.2.2 Lower Level 


The structural characteristics of the rotor are investigated at this level. The objectives are to 
minimize the tip displacements in cruise, in hover and in take-off. The most critical of these 
displacements are included in the formulation. In cruise, the elastic twist (<}> c ) and the inplane 
displacement (v c ) are critical. In hover and in take-off conditions, the vertical displacement (wh 
and w t> respectively) and the elastic twist (<J>h and 4> t , respectively) are significant. Therefore, these 
six displacements are selected as the individual objective functions to be minimized. Ply 
orientations are used as the discrete design variables at this level. Stress constraints are imposed 
and the Tsai-Wu failure criterion [22] is used to reduce the number of these constraints (see 
Chapter 4 for more details). 


5.3 Optimization Implementation 

The optimization strategy for this two level decomposition formulation is identical to the 
formulation used in Chapters 4. As before the Kreisselmeier-Steinhauser (K-S) function approach 
is used to formulate the multiple objective functions and the constraints into a single envelope 
function which is then extremized. Since only continuous design variables are used during 
optimization at the upper level, a nonlinear programming procedure (NLP) based on the Broyden- 
Fletcher-Goldfarb-Shanno (BFGS) algorithm is used. The lower level comprises only discrete 
design variables, therefore an optimization technique based on the simulated annealing algorithm is 
implemented at this level. Details of this formulation are found in Chapter 4. 

5.4 Results 

The same reference rotor is used in this study. For the hover and high speed cruise, the same 
operating conditions used in Chapters 3 and 4 are used. To simulate the take-off condition a load 
factor of 1.25 is used. Inclusion of the 12 percent down load effect, results in a take-off thrust of 
9100 lb. A rotational speed of 570 RPM is used and an altitude of 6695 feet is assumed to 
simulate a high altitude take-off. The flight conditions are summarized in Table 4. The blade is 
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discretized into 10 segments. The composite material used in the structural analysis is carbon- 
PEEK AS4/APC2 which has properties as listed in Table 2. 

At the upper level the design variables include the nodal values of the chord (c), twist (0), zero 
angle of attack (a z i) and thickness to chord ratio (t/c). The sweep distribution (A) is based on a 
quadratic lifting line. This yields a total of 45 design variables. The same scaling factor (k) used 
in Chapter 4 is employed in this objective to ensure that the box beam is maintained within the 
airfoil section (Eqn. 28). The tolerance on the lower level objective functions, used as optimal 
sensitivity parameters at this level, is initially set at 5 percent and is later relaxed to 20 percent for 
the tip bending displacements. An in-house code based on the Kreisselmeier-Steinhauser (K-S) 
function developed at Arizona State University is used as the optimization algorithm at this level. 
The search direction used during optimization is based on the BFGS algorithm and the two-point 
exponential expansion is used to approximate the objective functions and constraints. The K-S 
function multiplier 5 is initially set at 50 and increases during optimization to values as high as 210. 


Table 4 Summary of flight conditions 


Vehicle weight 

13,0001b 

Blade radius 

12.5 

Hover 

Altitude 

Sea level 

Thrust, Th 

7280 lb 

Rotational speed 

570 RPM 

Cruise 

Altitude 

25,000 ft 

Thrust, T c 

774 lb 

Rotational speed 

421 RPM 

Forward speed, Voo 

300 knots 

Take-off 

Altitude 

6695 ft 

Thrust, T t 
Rotational speed 

91001b 

570 RPM 


At the lower level, the design variables represent discrete values of the composite ply 
orientations. Since a symmetric and balanced lay-up is assumed in both the vertical and the 
horizontal walls, this leads to 12 independent design variables which can assume any one of the 7 
pre-selected values of ply angle orientations. A value of 5 is used for 8, the K-S function 
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multiplier. In the simulated annealing algorithm a value of 0.995 is used for the cooling rate, r c , 
and the initial temperature (To) is set to 1.0. 

In the multilevel problem, a total of 150 cycles is necessary for convergence at the upper level 
where a cycle consists of a “converged” design based upon one real analysis and several 
approximate function evaluations. A total of 5000 iterations are necessary at the lower level. Total 
convergence, including upper and lower levels, is achieved in five cycles. The results from this 
multiple design point optimization are presented in Tables 5 and 6 and Figs. 34-41. 

The upper level objective functions are presented in Table 5 and Fig. 34. The figures of merit 
in hover (FM h ) and in take-off (FM t ) are increased by 6.9 and 31 percent, respectively, from the 
reference values. A small increase (0.52%) is obtained in the cruise propulsive efficiency (Tic)- 
These trends can be explained by examining the rotor planform. 


Table 5 Summary of optimum results 



Reference 

Optimum 

— 

Objective functions 




Level 1 

FM (hover) 

0.760 

0.813 


T|ocx 

0.888 

0.893 


FM (take-off) 

0.617 

0.807 


Level 2 

wh (in) 

11.5 

6.95 


v c (in) 

-2.21 

-1.03 


w t (in) 

13.7 

7.91 


<t>h (deg) 

-2.05 

-1.77 


<|>c (deg) 

-1.22 

-1.19 

— 

<t>t (deg) 

- 1.73 

-1.39 


Solidity 

area weighted, O 

0.08075 

0.1055 

... 

thrust weighted, o t 

0.08913 

0.08976 
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A comparison of the optimum and reference chord distributions (c) are shown in Fig. 35 and 
significant differences must be noted. The optimum distribution closely resembles the well known 
ideal hover planform with notable exceptions at the root and at the tip. The proximity to the hover 
planform is due to the large solidity necessary to achieve the thrust required in both hover and in 
take-off. The deviation (from ideal hover planform) at the root is a result of an upper bound of 0.2 
which is imposed on the nondimensional chord (c/R) to avoid large chord sections. The deviation 
at the tip is due to a geometric constraint which is imposed to ensure that the box beam is 
maintained within the airfoil section. The reduction in the outboard section, relative to the 
reference rotor, is attributable to the lack of any maneuver margin requirement in the optimization 
problem formulation. Further, it is of interest to note that although the area-weighted solidity (a) is 
increased by nearly 3 1 percent, the thrust weighted solidity (Ot) is increased by only 7. 1 percent 
(Table 5). 

1 -, 

0.9- 
0 . 8 - 
0.7- 
0 . 6 - 
0.5- 
0.4- 
0.3- 
0 . 2 - 
0 . 1 - 
0 -- 

FMj, r| c FMt 

Figure 34 Summary of upper level objective functions 

The thickness to chord ratio (t/c) distributions of the reference and the optimum rotor 
configurations are presented in Fig. 36 and show large reductions from the reference rotor at the 
inboard sections of the optimum rotor. The thickness is also slightly reduced from the reference 


□ Reference 
| Multilevel 


0.52% 
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values at midspan locations and is slightly increased at the tip. The former is explained as an 
attempt to improve the rotor performance by reducing the profile drag. The latter represents the 
optimizer’s effort to satisfy the geometric constraint which ensures that the box beam is contained 
within the airfoil section. Since the optimum distribution is very similar to the reference 
distribution at midspan locations and slight increases are observed at the tip, the profile drag over 
the working section of the blade in case of high speed cruise is only slightly altered. This, coupled 
with the chord distribution over this section of the blade, results in only a slight improvement of 


the high speed cruise propulsive efficiency (r| c ). 


0.2 -n 
« 0.18- 
*0.16- 
1 0.14- 
« 0.12- 

1 01 ~ 

0 

0.08 - 
| 0.06 - 

1 0.04- 

Z 0.02 - 

0 -- 

0 
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Figure 35 Chord distributions 


The twist distributions (0) are presented in Fig. 37. The twist is reduced from reference to 
optimum values over the inboard section of the blade and is increased over the midspan section. At 
the outboard section of the blade, the twist remains almost unchanged after optimization. It must 
be noted from Fig. 37 that the largest differences between reference and optimum, in both inboard 
and midspan locations, are of the order of two degrees. However, the increase is achieved in a 
region of the blade which has greater resultant velocities. This reduces the collective pitch of the 
blade thereby reducing the overall angle of the attack of the blade. The result is a more even 
distribution of the angle of attack throughout the blade which subsequently reduces the drag 
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Thickness to chord ratio, t/c 


(Fig. 38a,b). The minimal changes in the angle of attack distribution in case of high speed cruise, 
between reference and optimum, partially explains the very small improvements obtained in r| c . 



0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 

Nondimensional radius, y/R 

Figure 36 Thickness to chord ratio distributions 

o 



Figure 37 Blade twist distributions 
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Figure 38a. Reference blade angle of attack Figure 38b. Optimum blade angle of attack 


Figure 39 shows the zero lift angle of attack distributions (a z i) for both the reference and the 
optimum rotors. Significant decreases are observed from the reference values except at the tip. 
The result is an increase in airfoil camber which improves the lift-to-drag ratio thereby improving 
performance. However, associated with increased camber is higher drag divergence Mach 
numbers (M<jd) which can adversely affect the cruise performance. Therefore to avoid large drag 
penalties caused by operating the blade at local Mach numbers above Mdd, the zero lift angle of 
attack is only slightly reduced at the tip. 

The rotor planform remains unswept after optimization. This is explained as follows. Since 
only static loading is included in this study and no aeroelastic stability requirements are imposed, 
the sweep of the blade was constrained during optimization such that the rotor could only be swept 
backwards. However, backward sweep induces large nose down pitching moments which 
increases the magnitude of the elastic twist which is used as an objective function and is minimized 
at each flight condition. As a result, the optimizer avoids sweeping the blades to remain in the 
feasible domain. 
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<t>h w h 

(hover) 


w h <j>c v c <|> t Wt 

,ver > (cruise) (take-off) 

Figure 40 Summary of lower level objective functions 


The lower level objective functions are presented in Table 5 and in Fig. 40 where large 
reductions from the reference rotor are observed in all six objective functions. The elastic twist in 
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hover OM, in cruise (W and in take-off <W are reduced by 14, 2.7 and 20 percent, respectively. 
The vertical displacement in hover <w„) is reduced by 40 percent and in take-off (w t ) by 42 
percent. The horizontal displacement in cruise (v c ) is reduced by 54 percent. 

The composite laminate stacking sequences are presented in Table 6. Note that in the reference 
blade, the horizontal and vertical walls are assumed to have the same stacking sequence. Further, 
since all of the laminates are considered to be symmetric about their midplane, only 12 of the 24 


total plies are presented in the table for each wall. The individual ply thickness used in this study ts 
0.001 in. which results in a total wall thickness of 0.24 in. The rearrangement of the stacking 
sequence after optimization represents a compromise between the conflicting requirements of 
reduced elastic twist and reduced transverse displacements. This is observed by noting the 
inclusion of ±30 degree plies in both the horizontal and the vertical walls. 

Table 6 Ply orientation angles 

Reference Optimum 

Horizontal Wall Vertical Wall 


Outer ply 


CP 
CP 
CP 
CP 
15° 
-15° 
15° 
-15° 
45° 
•45° 
45° 

midplane -45° 


15° 

-15° 

CP 

CP 

45° 

45° 

45° 

45° 

CP 

CP 

30P 

-30° 


15° 

-15° 

15° 

-15° 

30P 

-30° 

30P 

-30° 

CP 

CP 

15° 

-15° 


The elastic twist (*) distributions of the reference and optimum rotors are presented in Fig. 41. 
The figure shows that the elastic twist distributions are reduced in magnitude after optimization 
throughout the blade span in all three flight conditions. Note that the elastic twist values are all 
negative resulting in a nose down motion. Similar trends are observed in case of the transverse 
vertical displacement, the transverse horizontal displacement and the axial displacement. The 
reductions in all of the displacement distributions are attributable to the improved ply stacking 
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sequence and the increased thickness of the airfoil due to larger chords which in turn increases the 
box beam height. 



Radial location, y (ft) 
Figure 41 Elastic twist displacements 



A multiple design point optimization procedure was developed for the design of high speed 
prop-rotors. A multilevel decomposition optimization technique was used to decompose the 
problem into two levels. Aerodynamic performance was the objective of the upper level and the 
structural response was improved at the lower level. Optimization was performed simultaneously 
to include high speed cruise, hover and take-off. The Kreisselmeier-Steinhauser function was 
used to formulate the multiple objective optimization problems at each level. At the upper level, a 
nonlinear programming technique based on the Broyden-Fletcher-Goldfarb-Shanno method was 
used as the optimization algorithm. A simulated annealing algorithm was used for the discrete 
optimization problem at the lower level. A total of 5 global cycles were required for convergence. 
The optimum results were compared with a reference rotor. Following are some important 

observations. 
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1. The multilevel optimization procedure significantly improves the aerodynamic and 
structural response of the high speed prop-rotor blade at all three flight conditions. 

2. The chord distribution of the optimum rotor closely resembles the ideal hover chord 
distribution. Exceptions are noted at the root and at the tip where the optimizer is driven by 
the geometric constraints imposed on the problem. 

3 . The airfoil thickness to chord ratio is significantly decreased at inboard section of the blade 
after optimization. This reduces the profile drag over this portion of the blade which 
improves the performance. 

4. The optimum twist is reduced at inboard sections, increased over the mid section and 
remains largely unchanged at outboard sections. The result is a more optimum distribution 
of the local angle of attack which improves performance in hover and in take-off. The tip 
twist remains unchanged to maintain efficiency in high speed cruise. 

5 . The airfoil camber of the optimum blade is significantly increased from reference values to 
improve the airfoil lift-to-drag ratio resulting in increased aerodynamic efficiency. 

6. The simulated annealing algorithm successfully minimizes the tip displacements by altering 
the composite plate stacking sequences in the horizontal and vertical walls. The optimum 
composite stacking sequence represents a compromise between reduced elastic twist and 
reduced transverse deformation. This is manifested through the selection of ±30 degree 
plies in both the horizontal and the vertical walls. 

7 . A combination of improved stacking sequence and larger chord values leads to reduced 
elastic deformation in the optimum configurations. 


59 



6. Coupled Rotor/Wing Optimization 

Recently, research efforts have been initiated by Chattopadhyay et al. [4,5,15-18,29] to 
develop formal optimization techniques to address the conflicting issues involving tilt-rotor design. 
In Refs. 4 and 15 optimization procedures were developed to maximize the high speed cruise 
propulsive efficiency without degrading the hover figure of merit. The problem of individual blade 
aeroelastic stability in high speed cruise was included in the optimization formulation in Ref. 16. 
In Ref. 17, the drive system weight was minimized and the associated trade-off in cruise efficiency 
was investigated. The integrated aerodynamic, aeroelastic and structural optimization of the rotor 
was addressed in Ref. 5. In Ref. 18, a purely aerodynamic multiobjective optimization procedure 
was reported for improved high speed cruise and hovering performance using planform and airfoil 
characteristics as design variables. In Ref. 29, the aerodynamic and structural design criteria in 
both high speed cruise and hover were addressed by developing a multilevel decomposition based 
optimization procedure. In all of the above work, the rotor design was the primary concern and 
only the rotor planform parameters and geometric properties were included as design variables in 
the optimization formulation. 

The influence of the wing, however, is a critical design issue and requires in-depth study. 
Therefore, the development of a multidisciplinary optimization procedure involving simultaneous 
rotor-wing design is important in helping to understand the various instabilities which ultimately 
affect the aircraft performance. This study represents the results of a recent research effort aimed at 
studying the various trade-offs associated with such a design. In this study, the designs of the 
rotor and the wing are simultaneously addressed using formal optimization techniques. 

6. 1 Rotor/Wing Geometric Modeling 

The reference rotor is the representation of the XV- 15 rotor which was used in the previous 
objectives. The rotor is a three-bladed, gimballed rotor. The wing representation is based on the 
same aircraft configuration in which the rotor are mounted at the wing tips. The dynamic behavior 
of the wing is characterized by the first six frequencies, generalized masses and mode shapes of the 
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the load carrying member. The reference quantities correspond to the NASTRAN data of the XV- 
15 tilting prop-rotor aircraft [9,10,30]. 


6.1.1 Rotor Model 

The rotor planform characteristics are defined as follows. The chord, c( y ), twist angle of 
attack, 0(y) and the lifting line offset, x ac (y), are defined to have the following cubic spanwise 

distributions. 

0(y) = e L (y — 0.75) + 0 2 (4(y - 0.75) 2 - 3(y - 0.75)) 

+ 0 3 (15(y-O.75) 3 -2O(y-O.75) 2 +6(y-O.75)) (31) 

c(y) = c 0 + qy + c 2 (4y 2 - 3y) + c 3 (15y 3 - 20y 2 + 6y) (32) 

where 01-03 and C 0 -C 3 represent coefficients describing the spanwise variations of twist and chord, 
respectively. The offset in the twist distribution ( y-0.75) is used to ensure zero twist at 75 percent 
span. The shape functions used in these distributions are selected to maintain orthogonality 
between the functions over the range 0 to 1 which is beneficial to the optimization algorithm. 

A perturbational model, based on the chord distribution, is used to calculate the structural 
properties. The blade stiffness properties are assumed to vary from the reference values, which is 
the same baseline configurations as in the previous objectives, as follows. 


El 


xx 


El 


xx ref 


c(y) 

Cref(y) 


\ A 


(33) 


where EI XX is the flapping stiffness and c(y) is the chord. The lagging stiffness (EI ZZ ), the 
torsional rigidity (GJ) and the polar moment of inertia (Ie) are defined similarly. The fourth order 
relation is used since, in general, these quantities are based on fourth order dependence of the beam 
width and height. Assuming that these quantities are directly proportional to the chord, the fourth 
order dependence on the chord is obtained. The blade sectional mass and radius of gyration 
distributions, m( y) and My), respectively, are assumed to be quadratic as follows. 
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m 


( c(y) ^ 2 


m ref V c ref(y), 


and 




f c(y) ) 2 


L rref V 


Cref(y), 


(34) 


The quadratic dependence is assumed for these quantities since they tend to be dependent on the 
cross-sectional area, which again is assumed to directly proportional to the chord. Similarly, the 
center of gravity, tension center and lifting line offsets from the elastic axis are assumed to vary 
linearly with the chord ratio. 


6.1.2 Wing Model 

The aircraft wing is modeled as a cantilevered beam with tip masses representing the rotor 
system, transmission and engine pylon (Fig. 42). Although the reference wing planform is 
rectangular, the wing chord (C) is assumed to be linearly tapered to allow additional freedom to the 
optimizer as follows. 

C = c r [. + ;(*-!)] (35) 

where C r is the root chord, £ is the nondimensional wing span, X is taper ratio (X - C t /C r where C t 
is the tip chord). 

The load carrying structural member in the wing is modeled as a rectangular isotropic box beam 
with four independent wall thicknesses (Fig. 43). To reduce the number of design variables while 
keeping the design space sufficiently large, spanwise variations of wall thicknesses are assumed to 
be similar to the rotor chord and twist distributions as follows. 

t i(«;) = t io -+-t il <;-i-ti 2 (4<5 2 - 3<g) -t- tj 3 (15<g 3 -20<^ 2 36 

where the coefficients ti 0 - ti 3 describe the thickness distribution in each of the four walls 
0=1. 2, — , 4). 


6.2 Optimization Problem 

The coupled rotor/wing optimization problem is performed simultaneously under three different 
flight conditions. The first flight condition corresponds to sea level hover and the second flight 
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condition represents high speed cruise at an altitude of 25,000 feet. The third condition, which is 
used to investigate the aeroelastic stability of the coupled system, is a high speed windmill 
condition performed at an altitude of 25,000 feet. 

Z 




Figure 43 Wing cross section 
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6.2.1 Objective Functions 

The optimization problem addresses the multiple design requirements of maximizing the hover 
figure of merit (FM) and the propulsive efficiency in high speed cruise (r| c ) as well as minimizing 
the wing weight (W w j ng ) simultaneously. The wing weight is defined as follows, 
b 

W w ing = (1 + £) JpAdx (37) 

0 

where A is the cross-sectional area of the box beam, p is the density of the beam and b is the wing 
span. The factor e is used to include the effect of the nonstructural weight which is assumed to be 
one third of the structural weight. 

6.2.2 Constraints 

The following constraints are included in the optimization procedure to ensure efficient rotor 
and wing performance. To maintain rotor thrust at acceptable values during optimization, equality 
constraints are imposed on the thrust in hover and cruise. These constraints assume the following 
form. 

T h = T href and T c = T Cref (38) 

Additional constraints are imposed on the first natural frequency in hover (fj), the 
autorotational inertia (AI) and the blade weight (Wbiade)- To avoid ground resonance, the first 
natural frequency is constrained to be a prescribed amount above 1/rev. The autorotational inertia 
is constrained to be no less than 90 percent of the reference value to ensure that the optimum rotor 
can sufficiently autorotate in the event of engine failure. Since the autorotational inertia of the 
reference rotor is large, the minimum value of 90 percent of the reference value is selected to allow 
the optimizer more flexibility during the optimization process. An upper bound is imposed on the 
blade weight so that any reductions in the wing weight, which is an objective function, are not 
offset by an increase in the blade weight. These constraints are stated as follows. 

fl > (l+A)/rev (39) 
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AI > 0.9 AI re f 

Wblade — W re f 


(40) 

(41) 


where A = 0.2 represents a margin of safety which ensures that the first natural frequency in hover 
is above 1/rev. 

Since both the rotor and the wing designs are altered during optimization, it is important to 
impose aeroelastic stability constraints to prevent any destabilization of the rotor/wing in high 
speed cruise. Since the rotor is wing tip mounted for tilting prop-rotors, the first several modes of 
airframe, which are important for whirl flutter, are dominated by wing motion. To impose a factor 
of safety into the stability, the whirl flutter calculations are performed at the maximum dive speed 
which is 20 percent above the cruise condition. Also, the blade is trimmed to the windmill 
condition. These constraints are expressed as follows. 

A k < -x> k = 1, 2, , NMODE (42) 

where NMODE is the total number of modes considered in cruise and Ak is the real part of the 
stability root. The quantity \) denotes the minimum allowable blade damping and is defined to be 

equal to 0.01. 

Finally, constraints are imposed on the stresses at the root section of the wing box beam to 
ensure that the maximum stresses are below the allowable limit when multiplied by a factor of 
safety. These constraints assume the following form. 

ToctxFS < a a i (43) 

where NSEG is the number of wing segments, Toot is the octahedral stress at the wing root, FS is 
the factor of safety and a a i is the allowable stress. These stresses are calculated in both hover and 
cruise although they are more critical in hover. 

6.2.3 Design Variables 

To study the trade-offs associated with the simultaneous design of the rotor and the wing, 
design variables pertaining to both the rotor and the wing are used during optimization. The design 
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variables for the rotor include the coefficients which define the spanwise chord and twist 
distributions (co - C3 and 0 j - 63, respectively). As mentioned before, the shape functions used to 
describe the twist and the chord distributions are chosen to ensure orthogonality of the functions. 
Similarly the parameters that define the individual wall thicknesses in the wing-box (tjj) are used as 
design variables. The wing root chord (C r ) and taper ratio (A,) are also included as design 
variables. It must be noted that in order to ensure that the blade chord and wing thickness 
distributions are realistic (that is, positive throughout the span) it is necessary to further impose 
geometric constraints on these distributions. 

6.3 Analysis 

This section briefly describes the analysis procedures which are coupled within the 
optimization loop. This is followed by a description of the optimization technique. 

6.3.1 Rotor Analysis 

The analyses is performed using a comprehensive rotary wing analysis technique as 
implemented in the code CAMRAD [31]. The program calculates the section loading from the 
airfoil two-dimensional aerodynamic characteristics. It uses the lifting line or blade element 
approach and has corrections for yawed and three-dimensional flow effects. A wind tunnel trim 
option is used and the rotor is trimmed to a specific Cjtc value (Ct is the coefficient of thrust and 
a is the area-weighted solidity of the rotor) using the collective blade pitch. The value of Cjfa is 
calculated based on the prescribed values of thrust required at each flight condition. For the flutter 
calculations, a wind tunnel trim analysis is also performed. However, since this case is meant to 
be representative of a windmilling rotor, the rotor is trimmed to zero power. The aeroelastic 
stability analysis is performed using an assumption of axisymmetric flow. The blade is discretized 
into 51 structural elements and 24 aerodynamic elements. Six bending degrees of freedom, six 
torsional degrees of freedom and two gimbal degrees of freedom are used to model the rotor. The 
first six symmetric wing modes are also included in the whirl flutter calculations to yield a total of 
40 flutter modes. To reduce the computational effort, uniform inflow is assumed in all cases. The 
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blade response is calculated in CAMRAD using rotating free-vibration modes equivalent to a 
Galerkin analysis. The blade is trimmed at each optimization cycle so that an intermediate design, 
which is a feasible design, represents a trimmed configuration. 


6.3.2 Wing Analysis 

The wing aerodynamic and dynamic analysis is performed using an in-house developed code 
based on finite element analysis. The natural frequencies, generalized masses and generalized 
mode shapes of the wing are then used as inputs to describe the airframe, as required by 
CAMRAD. The load carrying member of the wing is modeled as a box beam with a rectangular — 
cross section (Fig. 43). Ten degrees of freedom are used to represent each beam element (Fig. 
42). To model the geometric offset between the rotor hub and the elastic axis of the wing box 
beam, rigid elements are included at the wing tip. The cross-sectional properties of every beam 
element are calculated using thin wall theory. To obtain a realistic wing representation, the 
reference configuration used in the optimization procedure is based on an untapered chord and wall — 
thickness distributions. 

The aerodynamic effects on the wing are evaluated using the lifting line theory with the 
assumption of quasi-steady aerodynamics which is used to compute the local lift coefficient 
distribution along the wing span. The aerodynamic forces are divided into a conservative part, 
which is analogous to a mechanical spring and a non conservative part which represents an _ 
aerodynamic damper. 

The derivation of the aerodynamic forces on a typical wing element is described here. 
Referring to Fig. 43 and assuming that the lift (L) acts at the quarter chord at a distance X a from 
elastic center, the lift can be expressed as follows. 


L = -ip_vic C|a 


a w +- 


w 


(44) 
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where is the air density, V M is the flight speed, C is wing chord length, q a is the lift curve 
slope, a w is the wing angle of attack and w is the velocity in the Z direction (Fig. 43). The work 
done by the aerodynamic force is then defined as follows. 


W = - [ L(w - a w X a )dX 

2 Jo 

where l e is the element length and w is the vertical displacement. Denoting 


(45) 


r = ^pooVlcc la 

and substituting into Eqn. 45 yields the following. 


(46) 


W = — 


if, 

2 Jo 


ra w (w X a a w )dX 


1 f r . 

2 Jo V “ 


w(w - X a a w )dX 


(47) 


The first integral in Eqn. 47 represents the conservative work, Wj and the second integral, W 2 , 
represents the non conservative work. Using a finite element approximation these quantities can 

be expressed as follows. 


Wj = — — x T K w x, W 2 — 2 * C w x 


(48) 


where x is the nodal displacement vector, K„ and C„ are analogous to the stiffness and damping 
matrices, respectively, normally associated with structural formulations and x is the time derivative 
of the nodal displacements. 


6.4 Optimization 

In this study, multiple objective functions are chosen and therefore the Kreisselmeier- 
Steinhauser (K-S) function is used to formulation the multiobjective optimization problem. Since 
the objective functions and the constraints must be evaluated several times before convergence is 
achieved, calculation of these values using exact analyses at each iteration is computationally 
prohibitive. Therefore the objective functions and constraints are approximated using a two-point 
exponential hybrid approximation technique. Details of these techniques are found in Chapters 

3-5. 
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6.5 Results 

The same reference rotor is again used for this study. The aerodynamic optimization is 
performed at a cruise altitude of 25,000 feet and a forward velocity of 300 knots with a rotational 
speed of 421 RPM. A vehicle weight of 13,000 pounds and, in this objective, an aircraft lift-to- 
drag ratio (L/D) of 5.3 are assumed. Therefore, the thrust in cruise is constrained to be 1226 
pounds for the two engine aircraft. In hover, the aircraft is assumed to be operating at sea level 
conditions with a rotational speed of 570 RPM. A 12 percent download effect from the rotor/wing 
interaction is used so that the thrust in hover is constrained to be 7280 pounds. The flutter 
calculations are performed at a speed 20 percent above the cruise speed to represent the maximum 
dive speed (360 knots) and an altitude of 25,000 feet with a rotational velocity of 421 RPM. The 
blade is trimmed to the windmill condition for the calculation of the flutter roots. The operating 
conditions are summarized in Table 7. 


Table 7 Summary of flight conditions 


Vehicle weight 

13,000 lb 

Blade radius 

12.5 

Hover 


Altitude 

Sea level 

Thrust, Th 

7280 lb 

Rotational speed 

570 RPM 

Cruise 


Altitude 

25,000 ft 

Thrust, T c 

1226 lb 

Rotational speed 

421 RPM 

Forward speed, 

300 knots 

Flutter 


Altitude 

25,000 ft 

Power, Cp 

0.0 

Rotational speed 

421 RPM 

Forward speed, V*, 

360 knots 


The results from the optimization are presented in Table 8 and Figs. 44 - 55. From Table 8 and 
Fig. 44, it is seen that the aerodynamic objective functions are both improved. The hover figure of 
merit is increased by 2.3 percent while the propulsive efficiency in high speed cruise is slightly 
increased (0.3 percent) after optimization. A comparison of the weights is presented in Table 8 and 
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Fig. 45 where it is seen that the wing weight (W w j ng ), which is also an objective function, is 
reduced by 14 percent from the reference value. The blade weight (Wbiade) which is included only 
as a constraint is also reduced significantly (6.3 percent) from the reference value. Further, the 
total rotor/wing weight which is defined for the three bladed rotor as 

W total = W wing + 3 W blade (49) 

is reduced by over 10 percent. These trends can be explained by examining the rotor and wing 
planforms. 


Table 8 Summary of optimum results 



Bounds 

Lower Upper 

Reference 

Optimum 

Objective functions 

FM 


0.7246 

0.7411 

Tic 


0.8475 

0.8502 

W w ing 


704.3 

606.1 

Constraints 

fi (per rev) 

1.200 

1.291 

1.327 

AI (lb-ft 2 ) 

3430 

3811 

3465 

W b iade Ob) 

213.4 

213.4 

200.0 

W t otal Ob) 

- 

1345 

1206 


The rotor blade chord distributions are presented in Fig. 46 and show that after optimization the 
blade chord has reduced values near the root and at the tip. This is due to the conflicting 
requirements posed by the hover and the cruise flight conditions. In hover, the rotor blade 
encounters negative angles of attack (a b ) at the tip which degrades performance (Fig. 47). The 
optimizer cannot simply reduce the twist at this location as this will degrade the high speed cruise 
propulsive efficiency. Similarly, the angle of attack distribution of the blade is negative near the 
root for the cruise condition (Fig. 48) and a simple increase in the twist at this location, to alleviate 
this problem, would degrade the hover figure of merit. The optimizer finds a compromise by 
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slightly reducing the twist near the root (Fig. 49) which favorably alters the collective trim angle, 
but more importantly by increasing the chord over the midspan region of the blade which is 
beneficial to both flight conditions. The chord is then reduced near the root to improve the high 
speed cruise propulsive efficiency and is reduced at the tip to improve the hover figure of merit. 



Figure 44 Optimum aerodynamic results 



Figure 45 Optimum structural results 
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Nondimcnsional chord, c/R 



Figure 46 Blade chord distributions 



Figure 47 Hover angle of attack distribution 
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Figure 48 High speed cruise angle of attack distribution 



Figure 49 Blade twist distributions 

The wing chord distributions are presented in Fig. 50. It is seen from this figure that the root 
chord remains unchanged and a very small amount of taper is introduced (A. = 0.99). Since the 

critical stresses occur at the root, the optimizer allows slight taper in the wing chord distribution 
which in turn tapers the dimensions of the box beam without exceeding the allowable stresses. If 
too much taper is introduced to reduce the weight, the stiffness of the beam will be reduced and the 
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aeroelastic stability will be degraded. This is why only a small taper is permissible in the wing 
chord distribution. 



Figure 50 Wing chord distributions 


The horizontal wall thickness of the wing box beam are presented in Fig. 51a. From the figure 
it is interesting to note that although the reference distributions for the top and the bottom walls are 
initially different, after optimization they are almost identical. The trend is to increase the 
thicknesses at the percent span. Near the tip, the thicknesses are slightly above the value at 75 
percent span location but are still drastically reduced from the reference values (28.5 and 22.5 
percent, respectively for the top and bottom walls). This is explained as follows. Near the root, 
the thicknesses are increased to satisfy the stress constraints which are initially violated. Since 
only the root stresses are included as constraints, the optimizer is now able to reduce the 
distributions throughout the span in order to reduce the wing weight. To ensure that the aeroelastic 
stability constraints are not violated (Fig. 52), the thicknesses near the tip are held nearly constant 
(relative to the 75 percent span value). The very small increases in the tip values, from the 75 
percent span location, are due to the nature of the orthogonal functions used in the formulation of 
the thickness distributions. If more higher order orthogonal functions were included in the 
formulation, the tip values would very likely equal the values at 75 percent span. It must be noted 
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that the cause of the initial constraint violation of the reference wing is due to the assumption of the 
box beam model. In order to closely match the natural frequencies, generalized masses and 
generalized mode shapes of the actual structural assembly which consists of stringers, spars and 
webs, the thickness distribution of the rectangular box beam produces stress constraints which are 
slightly violated initially. The authors acknowledge that a more accurate representation of the wing 
structural elements is necessary for further studies. 



Figure 51a Horizontal wall thicknesses distributions 
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The vertical wall thickness distributions are presented in Fig. 51b and shows that the optimum 
distributions are reduced from reference values throughout the blade span, including at the root. It 
must be noted here that the reference distributions in this case are the same for both the right and 
the left walls. The increased stiffnesses required to satisfy the stress constraints were achieved 
through an increase in the root thickness of the horizontal walls. Therefore, the optimizer is able to 
reduce the vertical wall thicknesses throughout the blade span to reduce the wing weight. It is of 
interest to note that despite the fact that cubic distributions are assumed for the vertical wall 
thicknesses, the optimum distributions in both cases essentially have only slight linear taper. 
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Figure 52 Aeroelastic stability characteristic exponents 

The first three mode shapes for the wing are presented in Figs. 53 - 55. The first mode 
(Fig. 53) is primarily an out-of-plane (flapping) motion. The second mode (Fig. 54) represents a 
primarily in-plane (chordwise bending) motion and the third mode (Fig. 55) is a elastic twist 
dominated motion. The figures also show the primary coupling associated with each mode. In 
Fig. 53 it is seen that the flapping motion of the first mode is slightly reduced from reference to 
optimum configuration. The coupling due to the pitching motion, however, is actually slightly 
increased in this mode. From Fig. 54 it is seen that the modes shapes for the second mode, which 
is primarily a lagging mode, are virtually unchanged after optimization. The third mode, which 
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represents a pitching dominated mode, is presented in Fig. 55 where it is again seen that primary 
mode is slightly reduced from reference to optimum, but the coupling, which in this case is due to 
the flapping motion, is again increased. The slight increases in the coupling of the wing modes are 
explained as follows. Wing weight minimization, one of the objectives, is accomplished through 
reductions of the wall thicknesses in the wing box beam and through the introduction of a slight 
wing chord taper. This results in small decreases in the wing stiffnesses and therefore small 
increases in the coupling of the mode shapes. A significant reduction in stiffnesses can increase 
whirl flutter instabilities. However, as shown in Fig. 52, although the aeroelastic stability 
constraints are made slightly more critical after optimization, they still remain stable. This is an 
important phenomenon and points to the critical trade-offs associated with optimum design for 
minimum wing weight while maintaining aeroelastic stability. 



Nondimensional wing span, Y/b 
Figure 53 First wing mode 
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Nondimensional wing span, Y/b 
Figure 54 Second wing mode 
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Figure 55 Third wing mode 


v6 Concluding Remarks 


A multidisciplinary optimization procedure is developed to analyze the coupled rotor/wing 


design for prop-rotor aircraft. The optimization procedure couples aerodynamic performance, 
aeroelastic stability and structural design for both high speed cruise and hover flight conditions. 


Objective functions include the maximization of the hover figure of merit and the high speed cruise 
propulsive efficiency and minimization of the wing weight. Constraints are imposed on the first 
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natural frequency in hover, the autorotational inertia, the blade weight and the maximum wing root 
stresses To prevent whirl flutter, constraints are also imposed on the real part of the stability roots. 
In addition, geometric constraints are also imposed. The optimization procedure is performed 
using the Kreisselmeier-Steinhauser (K-S) function technique and the Broyden-Fletcher-Goldfarb- 
Shanno (BFGS) algorithm. To reduce computational effort, an approximate analysis procedure 
based on the two-point exponential expansion and variable move limits is used. The results 
obtained are compared with an existing reference tilt-rotor configuration. The following important 
observations are made from this study. 

1 . Significant improvements are obtained in hover figure of merit and wing weight with small 
improvements in high speed cruise propulsive efficiency. The combined rotor/wing weight 
is also reduced significantly. 

2. The blade chord distribution is increased over the midspan region to exploit portions of the 
blade which simultaneously improve both the hover and high speed cruise performance. 
The chord distribution near the root and the tip are reduced due to the conflicting 
requirements posed by hover and cruise conditions. 

3. Wing weight reductions are achieved through reductions in the wall thicknesses although 
the horizontal walls are increased at the root to satisfy the root stress constraints which are 
violated in the reference design. The wing chord is only slightly tapered to maintain 
aeroelastic stability requirements. 

4. Reduction in wing weight results in small reductions in stiffness in the optimum 
configuration. This leads to increased coupling in wing modes. However, stability is still 
maintained in the optimum configuration. 
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Appendix - Nomenclature 


AI 

C 

Cp 

Ct 

c 

C0-C3 

Cd 

ci 

c la 

di 

EI.tx 

EIzz 

F 

E k 

* 

F k 

Fks 

fl 

fk 

FM 

Sj 

GJ 

H 

Ie 

J 

kij 

k r 

L 

M 

Mdd 

Mtip 


Roman 

Simulated annealing design variable vector 

Autorotational inertia, lb-ft 2 

Wing chord, ft 

Coefficient of power 

Coefficient of thrust 

Blade chord, ft 

Blade chord distribution parameters 
Coefficient of drag 
Coefficient of lift 
Lift curve slope 

Design variables vector of pre-selected values 
Lagging stiffness, lb-ft 2 
Flapping stiffness, lb-ft 2 
Axial force, lb 

Vector of original objective functions 

Vector of reduced objective functions 

K-S objective function 

First natural frequency in hover, per rev 

K-S function constraint vector 

Hover figure of merit 

Original constraint vector 

Torsional rigidity, lb-ft 2 

Interpolation matrix 

Polar moment of inertia, slug-ft 2 

Propeller advance ratio (jtVoo/^R) 

Stiffness matrix, p.s.i. 

Nondimensional blade radius of gyration 
Wing lift, lb 
Mach number 

Drag divergence Mach number 
Tip Mach number 


A-l 


My 

M z 

m 

NC 

NDV 

NF 

NSEG 

P 

P 

Pace 

Pn 

Q 

Qij 

Qy 

Qz 

R 

r 

r c 

T 

T k 

t 

tio‘U3 

t/c 

u, v, w 

Uj 

Voo 

w 

W blade 
W total 
^ wing 
VV 

X, y, Z 

X a 

x a/c 

y 


Flapping moment, lb-in 

Lagging moment, lb-in 

Blade sectional mass, slug/ft 

Number of constraints 

Number of design variables 

Number of objective functions 

Number of blade segments 

Rotor power, hp 

Random number, 0 < p < 1 

Acceptance probability 

Two-point exponential expansion exponent 

Rotor Torque 

Off-axis stiffness matrix for each lamina, p.s.i. 

Transverse horizontal shear force, lb 
Transverse vertical shear force, lb 
Blade radius, ft 
Radial location, ft 
Simulated annealing cooling rate 
Rotor thrust, lb 

Simulated annealing “temperature” 

Wall thickness, in 

Wall thickness distribution parameters, ft 
Thickness to chord ratio 
Elastic displacements 
Swirl velocity, ft/s 

Forward velocity, ft/s 2 
Resultant free stream velocity, ft/s 2 
Rotor blade weight, lb 
Combined wing/rotor weight, lb 
Wing weight, lb 
Wing flapping velocity, ft/s 
Reference axes 
Wing elastic axis offset, ft 
Lifting line offset 

Nondimensional radial location, y/R 
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Greek 


otb 

Blade angle of attack, deg. 

ai 

Local inflow angle, deg. 

CXw 

Wing angle of attack, rad. 

a z i 

Zero lift angle, deg. 

6 

K-S function draw down factor 

e 

Blade strain 

O 

Design variable vector 

<t> 

Elastic twist 

Yxn* YxX 

Shear strains 

Tl, C 

Local blade element axes 

ric-Tlt 

Rotor propulsive efficiency 

K 

Box beam scaling factor 

A 

Blade sweep, degrees 

X 

Taper ratio 

Ak 

Real part of the kth stability root 

X\v 

Out of plane warping function 


Helicopter advance ratio (Voo/OR) 

e 

Blade twist, degrees 

01-03 

Blade twist distribution parameters 

075 

Blade collective pitch angle, deg 

Poo 

Air density, slug/ft 3 

Cl, <52 

Material axis normal stresses, p.s. 

Gal 

Allowable stress, p.s.i. 

X12 

Shear stress, p.s.i. 

^oct 

Octahedral wing stress, p.s.i. 

Di 

Inflow velocity, ft/s 

Q 

Rotor angular velocity, r.p.m. 
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Nondimensional wing span, y/b 


c 


c 

h 

max 

ref 

t 


Subscripts 

High speed cruise 
Hover 

Maximum allowable 
Reference, or baseline, value 
Take-off 
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